Apparatus and program for analyzing distributed constant in a transmission line

ABSTRACT

The invention is aimed to analyze the characteristics of a transmission line only by inputting the specific distributed parameters without mesh dividing the transmission line to be analyzed into minimal unit required for the analysis. An arithmetic operation section  12  solves a differential equation regarding electromagnetic interaction of each line segment based on the predetermined distributed parameters, including a distributed self-inductance L i , a distributed resistance R i , and a distributed capacitance c i  and a distributed conductance g i  in relation to a reference potential plane in the transmission line, and a distributed mutual inductance M ij , a distributed capacitance C ij , and a distributed conductance G ij  between the transmission line and another transmission line, which are input from an input section  11 , and calculates the characteristic data of the transmission line.

TECHNICAL FIELD

The present invention relates to an apparatus for analyzing the distributed constant in a transmission line, and more particularly to an apparatus for analyzing the distributed constant, which is capable of analyzing the characteristics of the transmission line only by inputting the specific distributed parameters without mesh dividing the transmission line to be analyzed into minimal units required for the analysis. This apparatus for analyzing the distributed constant is especially suitable for a transmission line included in a spiral layered type distributed constant filter. This invention is also concerned with a program installed into an arithmetical operation section of the apparatus for analyzing the distributed constant.

BACKGROUND ART

Generally, a variety of high frequency devices are operable based on an electromagnetic interaction phenomenon described by the Maxwell's equation. Accordingly, a high frequency electromagnetic analysis is usually employed to analyze them. One of prior arts is well known as a design tool called Ansoft HFSS (High-Frequency Structure Simulator) developed by Ansoft Corp. in the United States.

This design tool developed by Ansoft Corp. is designed to involves automatically acquire the device characteristics of microwave, milliwave or wireless devices by solving the Maxwell's equation that is an operation principle for the high frequency devices, and can be used to determine the characteristics of antenna, transmission line, filter and connector by creating the shape of device to be analyzed using a three dimensional CAD interface, and setting up the material property values and an input/output port.

However, most of prior arts including the design tool developed by Ansoft Corp. make use of a Finite Element Method (FEM), and have the following problem:

In the Finite Element Method, the object having actually a complex shape and property is mesh divided into minimal units having a simple shape and property. Then, the characteristics of each element are approximately represented using a mathematical equation, and the simple equations are combined and solved simultaneously to predict the overall characteristics.

Therefore, the structure of an object to be analyzed must be mesh divided into simple minimal units required for the analysis, the units must be input into the three dimensional CAD, and the material property values for individual meshes must be specified. The procedure is time and labor intensive. Besides, the arithmetical operation spent several minutes to several hours in a super-computer class, and several days to several months in a workstation class, for example.

Accordingly, it is an object of the present invention to provide an apparatus and method for analyzing the distributed constant a transmission line to analyze the characteristics of the transmission line only by inputting the specific distributed parameters without mesh dividing the transmission line to be analyzed into minimal units required for the analysis, especially suitable for the transmission line included in a spiral layered type distributed constant filter.

SUMMARY OF THE INVENTION

In order to accomplish the above object, the present invention of claim 1 provides an apparatus for analyzing the distributed constant in a transmission line, comprising an input section for inputting the distributed parameters for each line segment virtually and arbitrarily divided in the transmission line to be analyzed, an arithmetic operation section for calculating the characteristic data of the transmission line by performing a predetermined arithmetic operation based on the distributed parameters input from the input section, and an output section for outputting the characteristic data, wherein the distributed parameters input from the input section include a distributed self-inductance L_(i), a distributed resistance R_(i), and a distributed capacitance c_(i) and a distributed conductance g_(i) in relation to a reference potential plane in the transmission line, and a distributed mutual inductance M_(ij), a distributed capacitance C_(ij), and a distributed conductance G_(ij) between the transmission line and another transmission line, wherein the arithmetic operation section calculates the characteristic data in the transmission line by solving a differential equation regarding electromagnetic interaction of the line segment based on the predetermined distributed parameters.

In the invention of claim 2, the arithmetic operation section comprises the secondary simultaneous differential equations as represented in the following formula 1 regarding electromagnetic interaction of the line segment with the voltage V(θ) of each line segment in the transmission line as an unknown variable (in the formula 1, V(θ) is a column vector of n×1 as represented in the following formula 2, Z is an impedance matrix created based on the distributed parameters as represented in the following formula 3, Y is an admittance matrix created based on the distributed parameters as represented in the following formula 4), in which the secondary simultaneous differential equations as represented in the formula 1 are solved by obtaining the eigenvalues γ_(i) ² (i=1 to n:n is the virtually divided number of line segments) of a complex coefficient matrix −ZY with the impedance matrix Z and the admittance matrix Y, thereby acquiring a general solution (p_(i), q_(i) are arbitrary constants, and u_(i) (i=1 to n) are eigenvectors corresponding to the eigenvalues γ_(i) ²) of voltage and current as represented in the following formula 5 and the following formula 6a or 6b at a limited position on the line segment. $\begin{matrix} {{\frac{\mathbb{d}^{2}\quad}{\mathbb{d}\theta^{2}}{V(\theta)}} = {{- {ZY}} \cdot {V(\theta)}}} & \left( {{Formula}\quad 1} \right) \\ {{V(\theta)} = \begin{pmatrix} {v_{1}(\theta)} \\ {v_{2}(\theta)} \\ \vdots \\ {v_{n}(\theta)} \end{pmatrix}} & \left( {{Formula}\quad 2} \right) \\ {Z = \begin{pmatrix} \left( {R_{1} + {j\quad\omega\quad L_{1}}} \right) & {j\quad\omega\quad M_{12}} & \cdots & {j\quad\omega\quad M_{1n}} \\ {j\quad\omega\quad M_{21}} & \left( {R_{2} + {j\quad\omega\quad L_{2}}} \right) & \cdots & {j\quad\omega\quad M_{2n}} \\ \vdots & \vdots & ⋰ & \vdots \\ {j\quad\omega\quad M_{n1}} & {j\quad\omega\quad M_{n2}} & \cdots & \left( {R_{n} + {j\quad\omega\quad L_{n}}} \right) \end{pmatrix}} & \left( {{Formula}\quad 3} \right) \\ {{{Y = \begin{pmatrix} \xi_{1} & \left( {G_{12} + {j\quad\omega\quad C_{12}}} \right) & \cdots & \left( {G_{1n} + {j\quad\omega\quad C_{1n}}} \right) \\ \left( {G_{21} + {j\quad\omega\quad C_{21}}} \right) & \xi_{2} & \cdots & \left( {G_{2n} + {j\quad\omega\quad C_{2n}}} \right) \\ \vdots & \vdots & ⋰ & \vdots \\ \left( {G_{n1} + {j\quad\omega\quad C_{n1}}} \right) & \left( {G_{n2} + {j\quad\omega\quad C_{n2}}} \right) & \cdots & \xi_{n} \end{pmatrix}}{{wherein},}}\quad} & \left( {{Formula}\quad 4} \right) \\ {{\xi_{i} = {- \left\{ {\left( {g_{j} + {j\quad\omega\quad c_{j}}} \right) + {\sum\limits_{{k = 1},{k \neq j}}^{\quad n}\left( {G_{jk} + {j\quad\omega\quad C_{jk}}} \right)}} \right\}}}{{V(\theta)} = {\sum\limits_{i = 1}^{n}{\left( {{p_{i}{\cosh\left( {\gamma_{i}\theta} \right)}} + {q_{i}{\sinh\left( {\gamma_{i}\theta} \right)}}} \right)u_{i}}}}} & \left( {{Formula}\quad 5} \right) \\ \begin{matrix} {{I(\theta)} = {Y{\int_{n}{{V(\theta)}{\mathbb{d}\theta}}}}} \\ {{= {Y{\sum\limits_{i = 1}^{n}{\frac{1}{\gamma_{i}}\left\{ {{p_{i}{\sinh\left( {\gamma_{i}\theta} \right)}} + {q_{i}{\cosh\left( {\gamma_{i}\theta} \right)}}} \right\} u_{1}\quad{or}}}}},} \end{matrix} & \left( {{Formula}\quad 6a} \right) \\ {{I(\theta)} = {{- Z^{- 1}}{\sum\limits_{i = 1}^{n}{\gamma_{i}\left\{ {{p_{i}{\sinh\left( {\gamma_{i}\theta} \right)}} + {q_{i}{\cosh\left( {\gamma_{i}\theta} \right)}}} \right\}{\mathbb{d}\theta}}}}} & \left( {{Formula}\quad 6b} \right) \end{matrix}$

In the invention of claim 3, the arithmetic operation section further calculates arbitrary constants (unknown variables) p_(i), q_(i) (i=1 to n) included in the general solution of voltage and current as represented in the formula 5 and the formula 6a or 6b under the boundary conditions of the line segment that are input from the input section.

In the invention of claim 4, the arithmetic operation section solves the secondary simultaneous differential equations as represented in the following formula 7 regarding electromagnetic interaction of the line segment with the current I(θ) of each line segment in the transmission line as an unknown variable, instead of the secondary simultaneous differential equations as represented in the formula 1, thereby acquiring a general solution of voltage and current as represented in the following formula 8 and the following formula 9a or 9b at a limited position on the line segment. $\begin{matrix} {{\frac{\mathbb{d}^{2}}{\mathbb{d}\theta^{2}}{I(\theta)}} = {{- {YZ}} \cdot {I(\theta)}}} & \left( {{Formula}\quad 7} \right) \\ {{I(\theta)} = {\sum\limits_{i = 1}^{n}{\left( {{p_{i}{\cosh\left( {\gamma_{i}\theta} \right)}} + {q_{i}{\sinh\left( {\gamma_{i}\theta} \right)}}} \right)u_{i}}}} & \left( {{Formula}\quad 8} \right) \\ \begin{matrix} {{V(\theta)} = {{- Z}{\int_{n}{{I(\theta)}{\mathbb{d}\theta}}}}} \\ {{= {{- Z}{\sum\limits_{i = 1}^{n}{\frac{1}{\gamma_{i}}\left\{ {{p_{i}{\sinh\left( {\gamma_{i}\theta} \right)}} + {q_{i}{\cosh\left( {\gamma_{i}\theta} \right)}}} \right\} u_{1}\quad{or}}}}},} \end{matrix} & \left( {{Formula}\quad 9a} \right) \\ {{V(\theta)} = {{- Z^{- 1}}{\sum\limits_{i = 1}^{n}{\gamma_{i}\left\{ {{p_{i}{\sinh\left( {\gamma_{i}\theta} \right)}} + {q_{i}{\cosh\left( {\gamma_{i}\theta} \right)}}} \right\} u_{i}}}}} & \left( {{Formula}\quad 9b} \right) \end{matrix}$

In the invention of claim 5, the arithmetic operation section comprises the linearly simultaneous differential equations as represented in the following formula 10 regarding electromagnetic interaction of the line segment with the voltage V(θ) and current I(θ) of each line segment in the transmission line as unknown variables (in the formula 10, W(θ) is a matrix containing V(θ)_(n×1) and I(θ)_(n×1) as represented in the following formula 11, A is a matrix as represented in the following formula 12 containing the impedance matrix Z created based on the distributed parameters as represented in the formula 3 and the admittance matrix Y created based on the distributed parameters as represented in the formula 4), in which the linearly simultaneous differential equations as represented in the formula 10 are solved by obtaining 2n eigenvalues γ_(i) ²(i=1 to 2n) of the matrix A, thereby acquiring a general solution (α_(i)(i=1 to 2n) are arbitrary constants, and u_(i) (i=1 to 2n) are eigenvectors corresponding to the eigenvalues γ_(i)) of voltage and current as represented in the following formula 13. $\begin{matrix} {{\frac{\mathbb{d}\quad}{\mathbb{d}\theta}{W(\theta)}} = {A \cdot {W(\theta)}}} & \left( {{Formula}\quad 10} \right) \\ {{W(\theta)} = {\begin{pmatrix} {V(\theta)}_{n \times 1} \\ {I(\theta)}_{n \times 1} \end{pmatrix} = \begin{pmatrix} {v_{1}(\theta)} \\ {v_{2}(\theta)} \\ \vdots \\ {v_{n}(\theta)} \\ {i_{1}(\theta)} \\ {i_{2}(\theta)} \\ \vdots \\ {i_{n}(\theta)} \end{pmatrix}}} & \left( {{Formula}\quad 11} \right) \\ {A = \begin{pmatrix} O_{n \times n} & {- Z} \\ Y & O_{n \times n} \end{pmatrix}_{2n \times 2n}} & \left( {{Formula}\quad 12} \right) \\ {{W(\theta)} = {\sum\limits_{i = 1}^{2n}{a_{i}{\mathbb{e}}^{{\gamma\mathbb{i}}\quad\theta}u_{i}}}} & \left( {{Formula}\quad 13} \right) \end{matrix}$

In the invention of claim 6, the arithmetic operation section sets up 2n simultaneous equations under the boundary conditions of the line segment input from the input section for each of the j-th (=1 to n) components of a voltage vector as represented in the following formula 14 and a current vector as represented in the following formula 15 in the transmission line, which are derived from the general solution as represented in the formula 13, thereby acquiring 2n arbitrary constants (unknown variables) α_(i). $\begin{matrix} {{v_{j}(\theta)} = {\sum\limits_{i = 1}^{2n}{\alpha_{i}{\mathbb{e}}^{{\gamma\mathbb{i}}\quad\theta}u_{ji}}}} & \left( {{Formula}\quad 14} \right) \\ {{i_{j}(\theta)} = {\sum\limits_{i = 1}^{2n}{\alpha_{i}{\mathbb{e}}^{{\gamma\mathbb{i}}\quad\theta}u_{n + {ji}}}}} & \left( {{Formula}\quad 15} \right) \end{matrix}$

In the invention of claim 7, the arithmetic operation section calculates an input/output voltage gain A_(v) from an input voltage and a terminating voltage of the signal line after obtaining the arbitrary constants of the general solution regarding the voltage and current.

In the invention of claim 8, the arithmetic operation section calculates an input impedance Z_(in) of the signal line based on a signal source voltage V_(s) and a signal source output impedance Z_(s) included in the boundary conditions given from the input section after obtaining the arbitrary constants of the general solution regarding the voltage and current.

In the invention of claim 9, the arithmetic operation section calculates an input/output voltage gain A_(v) from an input voltage and a terminating voltage of the signal line after obtaining the arbitrary constants of the general solution regarding the voltage and current, and calculates an input impedance Z_(in) of the signal line based on a signal source voltage V_(s) and a signal source output impedance Z_(s) included in the boundary conditions given from the input section, thereby acquiring an S parameter S₁₁ from the signal source output impedance Z_(s) and the input impedance Z_(in), and an S parameter S₂₁ from the S parameter S₁₁ and the input/output voltage gain A_(v).

Also, in the invention of claim 10, the signal line is involved in the equipment such as a connector, a sensor or a filter to allow the distributed constant characteristics of the equipment itself to be analyzed.

As defined in claim 11, the secondary simultaneous differential equations with the voltage V(θ) of each line segment as an unknown variable as represented in the formula 1 in claim 2 may be solved under the boundary conditions regarding the voltage of the line segment given from the input section in accordance with a numerical solution algorithm of the differential equation.

Similarly, as defined in claim 12, the secondary simultaneous differential equations with the current I(θ) of each line segment as an unknown variable as represented in the formula 7 in claim 4 may be solved under the boundary conditions regarding the current of the line segment given from the input section in accordance with a numerical solution algorithm of the differential equation.

Similarly, as defined in claim 13, the linearly simultaneous differential equations with the voltage V(θ) and current I(θ) of each line segment as unknown variables as represented in the formula 10 in claim 5 may be solved under the boundary conditions regarding the voltage and current of the line segment given from the input section in accordance with a numerical solution algorithm of the differential equation.

In this invention as defined in claim 14, the numerical solution algorithm is preferably any one of the Runge-Kutta methods, Shooting methods (Gun methods), and Relaxation methods.

Also, this invention comprises a program according to claim 15 for executing claim 2, a program according to claim 16 for executing claim 3, a program according to claim 17 for executing claim 4, a program according to claim 18 for executing claim 5, and a program according to claim 19 for executing claim 6.

Each of those programs is installed into the arithmetic operation section consisting of a computer, although its medium is not specifically limited. The program may be provided in a recording medium such as FD or CD-ROM, or downloaded via a communication line such as the Internet.

According to this invention, it is possible to analyze the transmission line characteristics only by inputting at least the distributed self-inductance L_(i), the distributed resistance R_(i), and the distributed capacitance c_(i) and the distributed conductance g_(i) in relation to the reference potential plane in the transmission line, and the distributed mutual inductance M_(ij), the distributed capacitance C_(ij), and the distributed conductance G_(ij) between the transmission line and another transmission line, as the distributed parameters for the transmission line to be analyzed, together with the predetermined boundary conditions of the transmission line, without mesh dividing the transmission line to be analyzed as in the Finite Element Method.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a block diagram showing a fundamental configuration of an apparatus for analyzing the distributed constant in a transmission line according to the present invention;

FIG. 2 is an explanatory view for defining a line segment in a relation with each pattern of the transmission line to be analyzed;

FIG. 3 is an equivalent circuit showing the distributed parameters that exist between two lines;

FIG. 4A is an exploded perspective view schematically showing a spiral layered type distributed constant device (filter) as analysis object;

FIG. 4B is an enlarged cross-sectional view of the essence of the spiral layered type distributed constant device;

FIG. 5 is an equivalent circuit diagram of the spiral layered type distributed constant device;

FIG. 6 is a matrix diagram showing a magnetic coupling state between lines in the spiral layered type distributed constant device;

FIG. 7 is a matrix diagram showing a capacitive coupling state in the spiral layered type distributed constant device;

FIG. 8 is a matrix diagram showing a magnetic coupling state between lines in the spiral layered type distributed constant device, when the length of one line segment is different;

FIG. 9 is a matrix diagram showing a capacitive coupling state between lines in the spiral layered type distributed constant device, when the length of one line segment is different;

FIG. 10 is a model circuit diagram that is simplified from the equivalent circuit diagram of FIG. 5;

FIG. 11 is a matrix diagram showing the simultaneous equations for deciding the unknown variables included in a general solution of a differential equation;

FIG. 12 is an equivalent circuit diagram of parallel three lines showing a specific analysis example of the invention;

FIG. 13A is a schematic diagram showing an output example of voltage distribution in the above analysis example;

FIG. 13B is a schematic diagram showing an output example of phase distribution in the above analysis example;

FIG. 14A is a schematic diagram showing another output example of phase distribution in the above analysis example;

FIG. 14B is a schematic diagram showing another output example of phase distribution in the above analysis example;

FIG. 15A is a graph showing an output example of the S₁₁ parameter in the above analysis example;

FIG. 15B is a graph showing an output example of the S₂₁ parameter in the above analysis example;

FIG. 15C is a graph showing the correlation in the size between S₁₁ parameter and S₂₁ parameter; and

FIG. 15D is a vector diagram showing the correlation in the phase between S₁₁ parameter and S₂₁ parameter.

DETAILED DESCRIPTION

First of all, an apparatus 10 for analyzing the distributed constant in a transmission line has a fundamental configuration of comprising an input section 11, an arithmetic operation section 12 and an output section 13, as shown in FIG. 1. The input section 11 may be a keyboard. The distributed parameters of the transmission line to be analyzed, including a distributed self-inductance L_(i)(H), a distributed resistance R_(i)(Ω), and a distributed capacitance c_(i)(F) and a distributed conductance g_(i)(S) in relation to a reference potential plane in the transmission line, and a distributed mutual inductance M_(ij)(H), a distributed capacitance C_(ij)(F), and a distributed conductance G_(ij)(S) between the transmission line and another transmission line, are input from the input section 11. Each of the distributed parameters is a complex number, except for the distributed resistance and the distributed conductance. The unit is in a unit system in which the line length is defined.

The arithmetic operation section 12 has a ROM 12 a and a working RAM 12 b, and the ROM 12 a stores a table for creating the differential equation or determinant. The arithmetic ability may be equal to or higher than that of the class of a personal computer commercially available on the market. The output section 13 may be preferably a color display or printer, and is connected to an external storage unit such as FD or CD-ROM, as needed.

This invention is applicable to most of the transmission lines, but in the case of a spiral transmission line as shown in FIG. 2, each turn has a different length, and a different line position for electromagnetic coupling, unlike a linear parallel distributed constant line.

Therefore, variable θ is used as the parameter. A variable range of θ is arbitrarily set. However, in a spiral layered type distributed constant device (filter) having the transmission line wound as shown in FIGS. 4A and 4B, due to the reasons that the line is electromagnetically coupled with another line for each turn and the distributed parameters are difficult to display in a non-dimensional number, it is assumed that the start point for each turn is 0 [rad] and the end point is 2π [rad] in the following.

When the start and end points are defined as above, the relation of coupling position is correctly kept by taking θ at right angles to the line in any shape of a parallel coupled line, a spiral line, a rectangular line and a heteromorphic loop line, as shown in FIG. 2. Accordingly, in the case of a concentric circle (approximately spiral), θ is the angle of the start point as seen from the center.

FIG. 3 shows a coupling distributed constant line including two lines 1 and 2. Referring to FIG. 3, a differential equation of distributed constant line will be described below. A voltage drop Δv_(i) with Δθ on line 1 is caused by a resistance R_(i) and a component of a self-inductance L₁ on the line itself and an induced voltage on the line 2, and represented by the following formula 100. Δν₁=−(R ₁ +jωL ₁)Δθi₁ −jωM ₁₂ Δθi ₂  (Formula 100)

Both sides of the formula 100 are divided by Δθ. Then, if Δθ goes to 0, the following formula 101 results. $\begin{matrix} {\frac{\mathbb{d}v_{1}}{\mathbb{d}\theta} = {{{- \left( {R_{1} + {j\quad\omega\quad L_{1}}} \right)}i_{1}} - {j\quad\omega\quad M_{12}i_{2}}}} & \left( {{Formula}\quad 101} \right) \end{matrix}$

The following formula 102 holds for the line 2 due to symmetry. $\begin{matrix} {\frac{\mathbb{d}v_{2}}{\mathbb{d}\theta} = {{{- j}\quad\omega\quad M_{21}i_{1}} - {\left( {R_{2} + {j\quad\omega\quad L_{2}}} \right)i_{2}}}} & \left( {{Formula}\quad 102} \right) \end{matrix}$

On the analogy from this result, if the voltage (potential, if seen from the reference plane) is represented by a function of position θ alone for n coupling distributed constant lines, the following formula 103 holds. $\begin{matrix} {{\frac{\mathbb{d}\quad}{\mathbb{d}\theta}{V(\theta)}} = {{- Z} \cdot {I(\theta)}}} & \left( {{Formula}\quad 103} \right) \end{matrix}$

Where V(θ) in the formula 103 is represented by the following formula 104, I(θ) is represented by the following formula 105, and Z is an impedance matrix that is a square matrix (here, M_(ij)=M_(ji)) represented by the following formula 106. $\begin{matrix} {{V(\theta)} = \begin{pmatrix} {v_{1}(\theta)} \\ {v_{2}(\theta)} \\ \vdots \\ {v_{n}(\theta)} \end{pmatrix}} & \left( {{Formula}\quad 104} \right) \\ {{I(\theta)} = \begin{pmatrix} {i_{1}(\theta)} \\ {i_{2}(\theta)} \\ \vdots \\ {i_{n}(\theta)} \end{pmatrix}} & \left( {{Formula}\quad 105} \right) \\ {Z = \begin{pmatrix} \left( {R_{1} + {j\quad\omega\quad L_{1}}} \right) & {j\quad\omega\quad M_{12}} & \cdots & {j\quad\omega\quad M_{1n}} \\ {j\quad\omega\quad M_{21}} & \left( {R_{2} + {j\quad\omega\quad L_{2}}} \right) & \cdots & {j\quad\omega\quad M_{2n}} \\ \vdots & \vdots & \ddots & \vdots \\ {j\quad\omega\quad M_{n1}} & {j\quad\omega\quad M_{n2}} & \cdots & \left( {R_{n} + {j\quad\omega\quad L_{n}}} \right) \end{pmatrix}} & \left( {{Formula}\quad 106} \right) \end{matrix}$

Then, a differential equation in view of the capacity between lines is derived. Firstly, the following formula 107 holds for i, like the differential equation regarding the voltage. $\begin{matrix} \begin{matrix} {{{- \Delta}\quad{i_{1}(\theta)}} = {\left( {G_{1} + {j\quad\omega\quad C_{1}}} \right)\Delta\quad{\theta \cdot \quad{v_{1}(\theta)}}}} \\ {{+ \left( {G_{12} + {j\quad\omega\quad C_{12}}} \right)}\Delta\quad{\theta \cdot \left( \quad{{v_{1}(\theta)} - {v_{2}(\theta)}} \right)}} \\ {{+ \left( {G_{13} + {j\quad\omega\quad C_{13}}} \right)}\Delta\quad{\theta \cdot \left( \quad{{v_{1}(\theta)} - {v_{3}(\theta)}} \right)}} \\ {\vdots} \\ {{+ \left( {G_{1n} + {j\quad\omega\quad C_{1n}}} \right)}\Delta\quad{\theta \cdot \left( \quad{{v_{1}(\theta)} - {v_{n}(\theta)}} \right)}} \end{matrix} & \left( {{Formula}\quad 107} \right) \end{matrix}$

In this formula 107, the first term is a leak current for the GND (ground), and the second term is a leak current due to a potential difference between line 1 and line 2 (hereinafter referred to as a leak current between adjacent lines). Where G₁₂, G₁₃, . . . G_(ln) are conductance [S/rad] between lines, and C₁₂, C₁₃, . . . C_(ln) are capacity [F/rad] between lines (G_(ij)=G_(ji), C_(ij)=C_(ji)). If this equation is represented in terms of the voltage for each line, the following formula 108 results. $\begin{matrix} {{{- \Delta}\quad{i_{1}(\theta)}} = {{\left\{ {\left( {G_{1} + {j\quad\omega\quad C_{1}}} \right) + {\sum\limits_{k = 2}^{n}\left( {G_{1k} + {j\quad\omega\quad C_{1k}}} \right)}} \right\}{{\Delta\theta} \cdot {v_{1}(\theta)}}} - {\sum\limits_{k = 2}^{n}{\left( {G_{1k} + {j\quad\omega\quad C_{1k}}} \right){{\Delta\theta} \cdot {v_{k}(\theta)}}}}}} & \left( {{Formula}\quad 108} \right) \end{matrix}$

Accordingly, both sides of the formula 108 are divided by −Δθ, and if Δθ goes to 0, the equation for each transmission line is generalized in the form of simultaneous equations, whereby the following formula 109 results. Where Y in the formula is an admittance matrix represented by the following formula 110 and ξ_(j) in the matrix Y is represented by the following formula 111. $\begin{matrix} {{\frac{\mathbb{d}\quad}{\mathbb{d}\theta}{I(\theta)}} = {Y \cdot {V(\theta)}}} & \left( {{Formula}\quad 109} \right) \\ {Y = \begin{pmatrix} \xi_{1} & \left( {G_{12} + {j\quad\omega\quad C_{12}}} \right) & \cdots & \left( {G_{1n} + {j\quad\omega\quad C_{1n}}} \right) \\ \left( {G_{21} + {j\quad\omega\quad C_{21}}} \right) & \xi_{2} & \cdots & \left( {G_{2n} + {j\quad\omega\quad C_{2n}}} \right) \\ \vdots & \vdots & \ddots & \vdots \\ \left( {G_{n1} + {j\quad\omega\quad C_{n1}}} \right) & \left( {G_{12} + {j\quad\omega\quad C_{n2}}} \right) & \cdots & \xi_{n} \end{pmatrix}} & \left( {{Formula}\quad 110} \right) \\ {\xi_{j} = {- \left\{ {\left( {g_{j} + {j\quad\omega\quad c_{j}}} \right) + {\sum\limits_{{k = 1},{k \neq j}}^{n}\left( {G_{jk} + {j\quad\omega\quad C_{jk}}} \right)}} \right\}}} & \left( {{Formula}\quad 111} \right) \end{matrix}$

The formula 103 is differentiated by θ and substituted for the formula 109, the secondary simultaneous differential equations for the parallel coupling distributed constant line represented by the following formula 112 (corresponding to formula 1 in claim 2) is obtained. $\begin{matrix} \begin{matrix} {{\frac{\mathbb{d}^{2}}{\mathbb{d}\theta^{2}}{V(\theta)}} = {{{- Z} \cdot \frac{\mathbb{d}\quad}{\mathbb{d}\theta}}{I(\theta)}}} \\ {= {{- {ZY}} \cdot {V(\theta)}}} \end{matrix} & \left( {{Formula}\quad 112} \right) \end{matrix}$

The secondary simultaneous differential equations are solved by obtaining n eigenvalues γ_(i) ² (i=1 to n) of a complex coefficient matrix −ZY. The result is represented by the following formula 113 (in this formula, u_(i) (i=1 to n) are eigenvectors corresponding to eigenvalues γ_(i) ²) $\begin{matrix} {{V(\theta)} = {{{\left( {{\alpha_{1}{\mathbb{e}}^{{+ \gamma_{1}}\theta}} + {\beta_{1}{\mathbb{e}}^{{- \gamma_{1}}\theta}}} \right){u_{1}\left( {{\alpha_{2}{\mathbb{e}}^{{+ \gamma_{2}}\theta}} + {\beta_{2}{\mathbb{e}}^{{- \gamma_{2}}\theta}}} \right)}u_{2}} + \ldots + {\left( {{\alpha_{n}{\mathbb{e}}^{{+ \gamma_{n}}\theta}} + {\beta_{n}{\mathbb{e}}^{{- \gamma_{n}}\theta}}} \right)u_{n}}} = {\sum\limits_{i = 1}^{n}{\left( {{\alpha_{i}{\mathbb{e}}^{{+ \gamma_{i}}\theta}} + {\beta_{i}{\mathbb{e}}^{{- \gamma_{i}}\theta}}} \right)u_{i}}}}} & \left( {{Formula}\quad 113} \right) \end{matrix}$

Moreover, substituting α_(i)+β_(i)=p_(i), α_(i)−β_(i)=q_(i) (p, q are arbitrary constants), this formula becomes an equation regarding the voltage using a hyperbolic function as represented by the following formula 114. $\begin{matrix} \begin{matrix} {{V(\theta)} = {\sum\limits_{i = 1}^{n}{\left( {{p_{i}\frac{{\mathbb{e}}^{{+ \gamma_{i}}\theta} + {\mathbb{e}}^{{- \gamma_{i}}\theta}}{2}} + {q_{i}\frac{{\mathbb{e}}^{{+ \gamma_{i}}\theta} - {\mathbb{e}}^{{- \gamma_{i}}\theta}}{2}}} \right)u_{i}}}} \\ {= \underset{\underset{n \times 1\quad{row}\quad{vector}}{︸}}{\sum\limits_{i = 1}^{n}{\left( {{p_{i}{\cosh\left( {\gamma_{i}\theta} \right)}} + {q_{i}{\sinh\left( {\gamma_{i}\theta} \right)}}} \right)u_{i}}}} \end{matrix} & \left( {{Formula}\quad 114} \right) \end{matrix}$

On the other hand, as an equation regarding the current, the following equation 115 is obtained from the formula 103. $\begin{matrix} \begin{matrix} {{I(\theta)} = {{- Z^{- 1}}\frac{\mathbb{d}\quad}{\mathbb{d}\theta}{V(\theta)}}} \\ {= {{- Z^{- 1}}\frac{\mathbb{d}\quad}{\mathbb{d}\theta}\left\{ {\sum\limits_{i = 1}^{n}{\left( {{p_{i}{\cosh\left( {\gamma_{i}\theta} \right)}} + {q_{i}{\sinh\left( {\gamma_{i}\theta} \right)}}} \right)u_{i}}} \right\}}} \\ {= {{- Z^{- 1}}{\sum\limits_{i = 1}^{n}{\gamma_{i}\left\{ {{p_{i}{\sinh\left( {\gamma_{i}\theta} \right)}} + {q_{i}{\cosh\left( {\gamma_{i}\theta} \right)}}} \right\} u_{i}}}}} \end{matrix} & \left( {{Formula}\quad 115} \right) \end{matrix}$

If the formula 109 regarding the current is indefinitely integrated, the following formula 116 equivalent to the formula 115 is obtained. Supposing that the current has no bias (DC) component, the integration constant is omitted. $\begin{matrix} \begin{matrix} {{I(\theta)} = {Y{\int{{V(\theta)}{\mathbb{d}\theta}}}}} \\ {= {Y{\int{\sum\limits_{i = 1}^{n}{\left( {{p_{i}{\cosh\left( {\gamma_{i}\theta} \right)}} + {q_{i}{\sinh\left( {\gamma_{i}\theta} \right)}}} \right)u_{i}{\mathbb{d}x}}}}}} \\ {= {Y{\sum\limits_{i = 1}^{n}{\left\{ {{p_{i}\frac{1}{\gamma_{i}}{\sinh\left( {\gamma_{i}\theta} \right)}} + {q_{i}\frac{1}{\gamma_{i}}{\cosh\left( {\gamma_{i}\theta} \right)}}} \right\} u_{i}}}}} \\ {= {Y{\sum\limits_{i = 1}^{n}{\frac{1}{\gamma_{i}}\left\{ {{p_{i}{\sinh\left( {\gamma_{i}\theta} \right)}} + {q_{i}{\cosh\left( {\gamma_{i}\theta} \right)}}} \right\} u_{i}}}}} \end{matrix} & \left( {{Formula}\quad 116} \right) \end{matrix}$

In this way, the secondary simultaneous differential equations of the formula 112 are solved to obtain a general solution for voltage and current as represented by the formulas 114 and 116. Accordingly, the j-th components of voltage and current vectors in the transmission line are represented by the following formulas 117 and 118. Using these formulas, 2n simultaneous equations are constructed under the boundary conditions to decide arbitrary constants (unknown variables) p_(i), q_(i). Where Ψ(i,j) in the formula 118 is given by the formula 119, and u_(ji) denotes the j-th element of the i-th eigenvector (n×1). $\begin{matrix} {{v_{j}(\theta)} = {{\sum\limits_{i = 1}^{n}{{\cosh\left( {\gamma_{i}\theta} \right)}{u_{ji} \cdot p_{i}}}} + {\sum\limits_{i = 1}^{n}{{\sinh\left( {\gamma_{i}\theta} \right)}{u_{ji} \cdot q_{i}}}}}} & \left( {{Formula}\quad 117} \right) \\ {{i_{j}(\theta)} = {{\sum\limits_{i = 1}^{n}{{\Psi\left( {i,j} \right)}{{\sinh\left( {\gamma_{i}\theta} \right)} \cdot p_{i}}}} + {\sum\limits_{i = 1}^{n}{{\Psi\left( {i,j} \right)}{{\cosh\left( {\gamma_{i}\theta} \right)} \cdot q_{i}}}}}} & \left( {{Formula}\quad 118} \right) \\ {{\Psi\left( {i,j} \right)} = {\frac{1}{\gamma_{i}}{\sum\limits_{k = 1}^{n}{y_{ik}u_{ki}}}}} & \left( {{Formula}\quad 119} \right) \end{matrix}$

In another form, the secondary simultaneous differential equations for the parallel coupling distributed constant line regarding the current may be constructed, using the vectors V(θ), I(θ) and the complex number matrixes Z and Y that are used in a secondary simultaneous differential equation model for the above formula 112 regarding the voltage (potential). These secondary simultaneous differential equations are shown in the following formula 120 (corresponding to the formula 7 in claim 4). $\begin{matrix} {{\frac{\mathbb{d}^{2}\quad}{\mathbb{d}\theta^{2}}{I(\theta)}} = {{- {YZ}} \cdot {I(\theta)}}} & \left( {{Formula}\quad 120} \right) \end{matrix}$

This equation, which has the same form as the formula 112, is solved by obtaining n eigenvalues γ_(i) ² (i=1 to n) for −YZ. The result is shown in the formula 121 (u_(i) (i=1 to n) are eigenvectors corresponding to the eigenvalues γ_(i) ², and p_(i), q_(i) are arbitrary constants). At this time, the voltage is represented by the formula 122. $\begin{matrix} {{I(\theta)} = {\sum\limits_{i = 1}^{n}{\left( {{p_{i}{\cosh\left( {\gamma_{i}\theta} \right)}} + {q_{i}{\sinh\left( {\gamma_{i}\theta} \right)}}} \right)u_{i}}}} & \left( {{Formula}\quad 121} \right) \\ \begin{matrix} {{V(\theta)} = {{- Z}{\int{{I(\theta)}{\mathbb{d}\theta}}}}} \\ {= {{- Z}{\sum\limits_{i = 1}^{n}{\frac{1}{\gamma_{i}}\left\{ {{p_{i}{\sinh\left( {\gamma_{i}\theta} \right)}} + {q_{i}{\cosh\left( {\gamma_{i}\theta} \right)}}} \right\} u_{i}}}}} \end{matrix} & \left( {{Formula}\quad 122} \right) \end{matrix}$

Accordingly, the j-th components of voltage and current vectors in the transmission line are represented by the following formulas 123 and 124. Where Ψ(i,j) in the formula 124 is given by the formula 125, and u_(ji) denotes the j-th element of the i-th eigenvector (n×1). $\begin{matrix} {{i_{j}(\theta)} = {{\sum\limits_{i = 1}^{n}{{\cosh\left( {\gamma_{i}\theta} \right)}{u_{ji} \cdot p_{i}}}} + {\sum\limits_{i = 1}^{n}{{\sinh\left( {\gamma_{i}\theta} \right)}{u_{ji} \cdot q_{i}}}}}} & \left( {{Formula}\quad 123} \right) \\ {{v_{j}(\theta)} = {{\sum\limits_{i = 1}^{n}{{\Psi\left( {i,j} \right)}{{\sinh\left( {\gamma_{i}\theta} \right)} \cdot p_{i}}}} + {\sum\limits_{i = 1}^{n}{{\Psi\left( {i,j} \right)}{{\cosh\left( {\gamma_{i}\theta} \right)} \cdot q_{i}}}}}} & \left( {{Formula}\quad 124} \right) \\ {{\Psi\left( {i,j} \right)} = {{- \frac{1}{\gamma_{i}}}{\sum\limits_{k = 1}^{n}{z_{jk}u_{ki}}}}} & \left( {{Formula}\quad 125} \right) \end{matrix}$

In still another form, the linearly simultaneous differential equations for the parallel coupling distributed constant line may be constructed, using the vectors V(θ), I(θ) and the complex number matrices ZY that are used in the secondary simultaneous differential equation model. The linearly simultaneous differential equations are shown in the following formula 126. $\begin{matrix} {{\frac{\mathbb{d}}{\mathbb{d}\theta}{W(\theta)}} = {A \cdot {W(\theta)}}} & \left( {{Formula}\quad 126} \right) \end{matrix}$

A vector W(θ) included in the formula 126 is represented by the following formula 127, and A is a matrix as represented by the following formula 128. $\begin{matrix} {{W(\theta)} = {\begin{pmatrix} {V(\theta)}_{n \times 1} \\ {I(\theta)}_{n \times 1} \end{pmatrix} = \begin{pmatrix} {\upsilon_{1}(\theta)} \\ {\upsilon_{2}(\theta)} \\ \vdots \\ {\upsilon_{n}(\theta)} \\ {i_{1}(\theta)} \\ {i_{2}(\theta)} \\ \vdots \\ {i_{n}(\theta)} \end{pmatrix}}} & \left( {{Formula}\quad 127} \right) \\ {A = \begin{pmatrix} 0_{n \times n} & {- Z} \\ Y & 0_{n \times n} \end{pmatrix}_{2n \times 2n}} & \left( {{Formula}\quad 128} \right) \end{matrix}$

The linearly simultaneous differential equation of the formula 126 is solved by obtaining n eigenvalues γ_(i) (i=1 to 2n) for the matrix A to obtain a general solution for the voltage and current as represented by the following formula 129 (u_(i) (i=1 to 2n) are eigenvectors corresponding to the eigenvalues γ_(i), and α_(i) (i=1 to 2n) are arbitrary constants). $\begin{matrix} {{W(\theta)} = {\sum\limits_{i = 1}^{2n}{\alpha_{i}{\mathbb{e}}^{\gamma\quad i\quad\theta}u_{i}}}} & \left( {{Formula}\quad 129} \right) \end{matrix}$

In this case, because 2n eigenvalues and eigenvectors must be obtained from 2n×2n matrix, the arithmetic operation amount is greater than solving the secondary simultaneous differential equations. However, the j-th components of voltage and current vectors in the transmission line are represented by the following formulas 130 and 131. Using these formulas, 2n simultaneous equations are constructed under the boundary conditions to decide 2n arbitrary constants (unknown variables) α_(i). $\begin{matrix} {{v_{j}(\theta)} = {\sum\limits_{i = 1}^{2n}{\alpha_{i}{\mathbb{e}}^{\gamma\quad i\quad\theta}u_{ji}}}} & \left( {{Formula}\quad 130} \right) \\ {{i_{j}(\theta)} = {{\sum\limits_{i = 1}^{2n}{\alpha_{i}{\mathbb{e}}^{\gamma\quad i\quad\theta}u_{n + {ji}}\quad j}} = {1\ldots\quad n}}} & \left( {{Formula}\quad 131} \right) \end{matrix}$

Then, the impedance matrix Z and the admittance matrix Y for use to solve the above differential equations will be described below, taking a distributed constant filter 20 as shown in FIG. 4 as an example. FIG. 4A is an exploded perspective view, and FIG. 4B is an enlarged cross-sectional view showing a layered state.

This distributed constant filter 20 is a spiral layered type distributed constant filter containing a filter main body 22 having the transmission line wound like a spiral within the pot cores 21, 21 of Mn-Zn ferrite, the filter main body 22 comprising two tape shaped dielectric lines 23 a, 23 b made of polyimide, for example.

On one side of a tape shaped dielectric line 23 a, a transmission line 24 composed of a copper foil having a thickness of 35 μm, for example, is formed. Also, on one side of a tape shaped dielectric line 23 b, a ground (GND) line 25 composed of a copper foil having a thickness of 35 μm is formed. These two tape shaped dielectric lines 23 a, 23 b are wound like spiral. FIG. 5 shows an equivalent circuit of the distributed constant filter 20.

Herein, it is assumed that each line 24, 25 is virtually divided for each turn into 2n line segments. That is, supposing that each line 24, 25 is linked continuously, the transmission line 24 has the line segments 1 to n allocated, and the grounding line 25 has the line segments n+1 to 2n allocated. The line segment 1 of the transmission line 24 and the line segment n+1 of the grounding line 25 are associated, and the line segment n of the transmission line 24 and the line segment 2n of the grounding line 25 are associated, as shown in FIG. 5.

As shown in FIG. 5, there are various distributed parameters in each line segment. Since the magnetic coupling exists in all the line segments (each turn), the coupling between the distributed self-inductance L_(i) [H/rad] and the distributed mutual inductance M_(ij) [H/rad] is represented by a matrix as shown in FIG. 6.

In the matrix of FIG. 6, the diagonal element L_(i) indicates the self-inductance in each line segment (turn), and the off-diagonal element M_(ij) indicates the mutual coupling between line i and line j. The matrix of each ¼ block represents a magnetic coupling state in the transmission line 24 in the left upper ¼ block, and a magnetic coupling state in the grounding line 24 in the right lower ¼ block. And the remaining matrixes of ¼ block represent the magnetic coupling state between the transmission line 24 and the grounding line 25 in the left lower and right upper blocks. This matrix is multiplied by jω, and the resistance components R, are added to the diagonal components, whereby the impedance matrix Z as represented by the formula 106 is constructed.

On the other hand, the capacitive coupling exists only between adjacent lines via the tape shaped dielectric line, and is represented by a matrix as shown in FIG. 7, assuming that the distributed capacity between lines is C_(ij) [F/rad], and the distributed capacity of each line in relation to the reference potential plane is c_(i) [F/rad]. The admittance matrix Y as represented by the formula 110 is created based on this matrix.

The division and allocation of line segments may not be made for each turn. For example, assuming that one line segment is 0.5 turn, when one line is divided into n, double the n value in the above case is attained, whereby the magnetic coupling is in a combination of (even number segment) and (even number segment), or (odd number segment) and (odd number segment). At this time, the magnetic coupling matrix is represented in FIG. 8, and the capacitive coupling matrix is represented as shown in FIG. 9.

In this way, in virtually dividing one line, according to the invention, the line may be divided for any number of turn. If the distributed parameters for the coupling relation between line segments are input from the input section 11, the arithmetic operation section 12 solves the differential equations by combining the above matrixes.

In another form of the invention, the secondary simultaneous differential equations with the voltage V(θ) of each line segment as an unknown variable, as represented in the formula 112 (formula 1), is solved in accordance with a numerical solution algorithm for the differential equation under the boundary conditions regarding the voltage of each line segment that are given from the input section 11.

Similarly, the secondary simultaneous differential equations with the current I(θ) of each line segment as an unknown variable, as represented in the formula 120 (formula 7), is solved in accordance with a numerical solution algorithm for the differential equation under the boundary conditions regarding the current of each line segment that are given from the input section 11.

Also, similarly, the linearly simultaneous differential equations with the voltage V(θ) and current I(θ) of each line segment as unknown variables, as represented in the formula 126 (formula 10), is solved in accordance with a numerical solution algorithm for the differential equation under the boundary conditions regarding the voltage and current of each line segment that are given from the input section 11.

The numerical solution algorithm is preferably any one of the Runge-Kutta methods, Shooting methods (Gun methods), and Relaxation methods. These numerical solution algorithms involve providing the boundary conditions (boundary values) at predetermined positions in the line segment θ=Ai (1≦i), and acquiring the value at a desired position θ=B to be consistent with the boundary values at these points. Thereby, the solution values can be approximately obtained directly from the differential equation.

Referring to FIG. 10, a way of acquiring arbitrary variables (unknown variables) p_(i), q_(i) contained in the general solution of the differential equation will be described below. FIG. 10 is an equivalent circuit diagram that is simplified from the equivalent circuit of FIG. 5. In deciding the unknown variables p_(i), q_(i), the boundary conditions as listed in Table 1 are set up by the input section 11 in this example. V_(s): Signal source voltage (arbitary) Z_(s): Signal source output impedance Z_(g1): Signal source side grounding impedance Z_(g2): Terminating side grounding impedance Z_(t): Terminating load impedance Z_(r): Transmission line/grounding line r-1 turn connection impedance Z_(gi·r): Transmission line/grounding line r-1 turn output side grounding impedance Z_(gi·r): Transmission line/grounding line r-1 turn input side grounding impedance Z_(li): Impedance with lead wire at input end of transmission line Z_(lo): Impedance with lead wire at ouput end of transmission line

The arithmetic operation section 12 sets up the following formulas 132 to 137 under the given boundary conditions in accordance with the Kirchhoff's laws. ν₁(0)+(Z _(a) +Z1 _(i))·i₁(0)=V _(a)  (Formula 132) i_(n+1)(0)+ν_(n+1)(0)/Z _(gl)=0  (Formula 133) ν_(n)(2π)−(Z _(t) +Z1 _(o))·i _(n)(2π)=0  (Formula 134) σ_(2n)(2π)−Z _(g2) ·i _(2n)(2π)=0  (Formula 135) i _(r)(0)−i _(r−1)(2π)+ν_(r)(0)/Z _(go·r)+ν_(r−1)(2π)/Z _(gi·r)=0(r=2 . . . n,n+2 . . . 2n)  (Formula 136) (ν_(r−1)(2π)−ν_(r)(0))/Z _(r) −i _(r)(0)−ν_(r)(0)/Z _(go·r)=0(r=2 . . . n,n+2 . . . 2n)  (Formula 137)

Among them, the formulas 132 and 133 are boundary condition expressions regarding the signal source V_(s), and the formulas 134 and 135 are boundary condition expressions regarding the terminating load Z_(t). Also, the formulas 136 and 137 are boundary condition expressions at a connection point of each line. In each of the formulas, r takes a value from 2 to n and from n+2 to 2n, and 4n equations are prepared from all the formulas 132 to 137 to acquire the 4n unknown constants.

In the device as shown in FIG. 10, V_(s), Z_(s) and Z_(t) are appropriately given according to the conditions, and supposing a structure in which the grounding line 25 is connected to the GND at the input side, it follows that Z_(gi)=0 and Z_(g2)=∞. Also, since the transmission line 24 is continuously linked, it is sufficient that Z_(r)=0 and Z_(gi·r)=Z_(g0·r)=∞. In the following, the operation of the arithmetic operation section (computer) 12 to obtain the unknown variables will be described.

First of all, the formulas 132 and 133 that are boundary condition expressions regarding the signal source V_(s) are performed by the following formulas 138 and 139. $\begin{matrix} {{{\upsilon_{1}(0)} + {Z_{s} \cdot {i_{1}(0)}} + {{Zl}_{i} \cdot {i_{1}(0)}}} = {{{\sum\limits_{i = 1}^{2n}{u_{1i} \cdot p_{i}}} + {\left( {Z_{s} + {Zl}_{i}} \right)\left( {\sum\limits_{i = 1}^{2n}{{\Psi\left( {i,1} \right)} \cdot q_{i}}} \right)}} = {{{\sum\limits_{i = 1}^{2n}{\underset{\underset{f_{{\upsilon sp}\quad{(i)}}}{︸}}{u_{1i}} \cdot p_{i}}} + {\sum\limits_{i = 1}^{2n}{\underset{\underset{f_{{\upsilon sq}{(i)}}}{︸}}{\left\{ {\left( {Z_{s} + {Zl}_{i}} \right) \cdot {\Psi\left( {i,1} \right)}} \right\}} \cdot q_{i}}}} = V_{s}}}} & \left( {{Formula}\quad 138} \right) \\ {{{{\upsilon_{n + 1}(0)}/Z_{g1}} + {i_{n + 1}(0)}} = {{{\sum\limits_{i = 1}^{2n}{{u_{n + {1 \cdot i}}/Z_{g1}} \cdot p_{i}}} + {\sum\limits_{i = 1}^{2n}{{\Psi\left( {i,{n + 1}} \right)} \cdot q_{i}}}} = {{{\sum\limits_{i = 1}^{2n}{\underset{\underset{f_{{isp}{(i)}}}{︸}}{u_{n + {1 \cdot i}}/Z_{g1}} \cdot p_{i}}} + {\sum\limits_{i = 1}^{2n}\underset{\underset{f_{{isq}{(i)}}}{︸}}{{\Psi\left( {i,{n + 1}} \right)} \cdot q_{i}}}} = 0}}} & \left( {{Formula}\quad 139} \right) \end{matrix}$

The formulas 134 and 135 that are boundary condition expressions regarding the terminating load Z_(t) are performed by the following formulas 140 and 141. $\begin{matrix} {{{\upsilon_{n}\left( {2\pi} \right)} - {\left( {Z_{t} + Z_{l_{o}}} \right) \cdot {i_{n}\left( {2\pi} \right)}}}{\begin{matrix} {\quad{= {{\sum\limits_{i = 1}^{2n}{{\cosh\left( {2\quad\pi\quad\gamma_{i}} \right)}{u_{ni} \cdot p_{i}}}} + {\sum\limits_{i = 1}^{2n}{\sin\quad{h\left( {2{\pi\gamma}_{i}} \right)}{u_{ni} \cdot q_{i}}}} - \left( {Z_{t} + {Zl}_{o}} \right)}}} \\ {\left( {{\sum\limits_{i = 1}^{2n}{{\Psi\left( {i,n} \right)}\sin\quad{{h\left( {2\pi\quad\gamma_{i}} \right)} \cdot p_{i}}}} + {\sum\limits_{i = 1}^{2n}{{\Psi\left( {i,n} \right)}\cos\quad{{h\left( {2{\pi\gamma}_{i}} \right)} \cdot q_{i}}}}} \right)} \end{matrix}\begin{matrix} {\quad{= {{\sum\limits_{i = 1}^{2n}{\underset{\underset{f_{gtp}{(i)}}{︸}}{\left( {{\cos\quad{h\left( {2{\pi\gamma}_{i}} \right)}u_{ni}} - {{\left( {Z_{t} + {Zl}_{o}} \right) \cdot {\Psi\left( {i,n} \right)}}\sin\quad{h\left( {2{\pi\gamma}_{i}} \right)}}} \right)} \cdot p_{i}}} +}}} \\ {= {\sum\limits_{i = 1}^{2n}{\underset{\underset{f_{gtq}{(i)}}{︸}}{\left( {{\sin\quad{h\left( {2{\pi\gamma}_{i}} \right)}u_{ni}} - {{\left( {Z_{t} + {Zl}_{o}} \right) \cdot {\Psi\left( {i,n} \right)}}\cos\quad{h\left( {2{\pi\gamma}_{i}} \right)}}} \right)} \cdot q_{i}}}} \\ {= 0} \end{matrix}}{{\upsilon_{2n}\left( {2\pi} \right)} - {Z_{g2} \cdot {i_{\quad_{2n}}\left( {2\pi} \right)}}}} & \left( {{Formula}\quad 140} \right) \\ \begin{matrix} {\quad{= {{\sum\limits_{i = 1}^{2n}{\cos\quad{h\left( {2\quad\pi\quad\gamma_{i}} \right)}{u_{2{n \cdot i}} \cdot p_{i}}}} + {\sum\limits_{i = 1}^{2n}{\sin\quad{h\left( {2{\pi\gamma}_{i}} \right)}{u_{2{n \cdot i}} \cdot q_{i}}}} -}}} \\ {Z_{2g}\left( {{\sum\limits_{i = 1}^{2n}{{\Psi\left( {i,{2n}} \right)}\sin\quad{{h\left( {2\pi\quad\gamma_{i}} \right)} \cdot p_{i}}}} + {\sum\limits_{i = 1}^{2n}{{\Psi\left( {i,{2n}} \right)}\cos\quad{{h\left( {2\pi\quad\gamma_{i}} \right)} \cdot}}}} \right.} \\ {{\left. q_{i} \right){\sum\limits_{i = 1}^{2n}{\underset{\underset{f_{gtp}{(i)}}{︸}}{\left( {{\cos\quad{h\left( {2{\pi\gamma}_{i}} \right)}u_{2{n \cdot i}}} - {{Z_{g2} \cdot {\Psi\left( {i,{2n}} \right)}}\sin\quad{h\left( {2{\pi\gamma}_{i}} \right)}}} \right)} \cdot p_{i}}}} +} \\ {\sum\limits_{i = 1}^{2n}{\underset{\underset{f_{gtq}{(i)}}{︸}}{\left( {{\sin\quad{h\left( {2{\pi\gamma}_{i}} \right)}u_{2{n \cdot i}}} - {{Z_{g2} \cdot {\Psi\left( {i,{2n}} \right)}}\cos\quad{h\left( {2{\pi\gamma}_{i}} \right)}}} \right)} \cdot q_{i}}} \\ {= 0} \end{matrix} & \left( {{Formula}\quad 141} \right) \end{matrix}$

Also, the formulas 136 and 137 that are boundary condition expressions at a connection point of each transmission line are performed by the following formulas 142 and 143. $\begin{matrix} {{{i_{r}(0)} - {i_{r - 1}\left( {2\pi} \right)} + {{\upsilon_{r}(0)}/Z_{{go} \cdot r}} + {{\upsilon_{r - 1}\left( {2\pi} \right)}/Z_{{gi} \cdot r}}}\begin{matrix} {\quad{= {{\sum\limits_{i = 1}^{2\pi}{{\Psi\left( {i,r} \right)} \cdot q_{i}}} - {\underset{i = 1}{\overset{2n}{\left( \sum \right.}}{\Psi\left( {i,{r - 1}} \right)}\sin\quad{{h\left( {2{\pi\gamma}_{i}} \right)} \cdot p_{i}}} + {\sum\limits_{i = 1}^{2n}{\Psi\left( {i,{r - 1}} \right)\cos\quad{{h\left( {2{\pi\gamma}_{i}} \right)} \cdot}}}}}} \\ {\left. q_{i} \right) + {\frac{1}{Z_{{go} \cdot r}}\left( {\sum\limits_{i = 1}^{2n}{u_{ri} \cdot p_{i}}} \right)} + {\frac{1}{Z_{{gi} \cdot r}}\left( {{\sum\limits_{i = 1}^{2n}{\cos\quad{h\left( {2{\pi\gamma}_{i}} \right)}{u_{r - {1 \cdot i}} \cdot p_{i}}}} +} \right.}} \\ \left. {\sum\limits_{i = 1}^{2n}{\sin\quad{h\left( {2{\pi\gamma}_{i}} \right)}{u_{r - {1 \cdot i}} \cdot q_{i}}}} \right) \\ {= {{\sum\limits_{i = 1}^{2n}{\underset{\underset{f_{\upsilon\quad{lq}}{({r,i})}}{︸}}{\left. \left( {{{- {\Psi\left( {i,{r - 1}} \right)}}\sin\quad{h\left( {2{\pi\gamma}_{i}} \right)}} + {\frac{1}{Z_{{go} \cdot r}}u_{ri}} + {\frac{1}{Z_{{gi} \cdot r}}\cos\quad{h\left( {2{\pi\gamma}_{i}} \right)}u_{r - {1 \cdot i}}}} \right) \right)} \cdot p_{i}}} +}} \\ {\sum\limits_{i = 1}^{2n}{\underset{\underset{f_{vlq}{({r,i})}}{︸}}{\left( {{\Psi\left( {i,r} \right)} - {{\Psi\left( {i,{r - 1}} \right)}\cos\quad{h\left( {2{\pi\gamma}_{i}} \right)}} + {\frac{1}{Z_{{gi} \cdot r}}\sin\quad h\quad\left( {2{\pi\gamma}_{i}} \right)u_{r - {1 \cdot i}}}} \right)} \cdot q_{i}}} \\ {= 0} \end{matrix}} & \left( {{Formula}\quad 142} \right) \\ {{{{- {\upsilon_{r}(0)}}/Z_{r}} + {{\upsilon_{r - 1}\left( {2\pi} \right)}/Z_{r}} - {i_{r}(0)} - {{\upsilon_{r}(0)}/Z_{{go} \cdot r}}}\begin{matrix} {\quad{= {{{- \frac{1}{Z_{r}}}\left( {\sum\limits_{i = 1}^{2n}{u_{ri} \cdot p_{i}}} \right)} + {\frac{1}{Z_{r}}\left( {{\sum\limits_{i = 1}^{2n}{\cos\quad{h\left( {2{\pi\gamma}_{i}} \right)}{u_{r - {1 \cdot i}} \cdot p_{i}}}} +} \right.}}}} \\ {\left. {\sum\limits_{i = 1}^{2n}{\sin\quad{h\left( {2{\pi\gamma}_{i}} \right)}{u_{r - {1 \cdot i}} \cdot q_{i}}}} \right) - {\sum\limits_{i = 1}^{2n}{{\Psi\left( {i,r} \right)} \cdot q_{i}}} - {\frac{1}{Z_{{go} \cdot r}}\left( {\sum\limits_{i = 1}^{2n}{u_{ri} \cdot p_{i}}} \right)}} \\ {= {{\sum\limits_{i = 1}^{2n}{\underset{\underset{f_{{ilp}{({r,i})}}}{︸}}{\left( {{{- \frac{1}{Z_{r}}} \cdot u_{ri}} + {\frac{1}{Z_{r}}\cos\quad{h\left( {2{\pi\gamma}_{i}} \right)}u_{r - {1 \cdot i}}} - {\frac{1}{Z_{{go} \cdot r}} \cdot u_{ri}}} \right)} \cdot p_{i}}} +}} \\ {\sum\limits_{i = 1}^{2n}{\underset{\underset{f_{{ilq}{({r,i})}}}{︸}}{\left( {{\frac{1}{Z_{r}}\sin\quad{h\left( {2{\pi\gamma}_{i}} \right)}u_{r - {1 \cdot i}}} - {\Psi\left( {i,r} \right)}} \right)} \cdot q_{i}}} \\ {= 0} \end{matrix}} & \left( {{Formula}\quad 143} \right) \end{matrix}$

Thereby, it follows that 4n equations are prepared for 4n unknown variables p₁ to p_(n), q₁ to q_(n), p_(n+1) to p_(2n), and q_(n+1) to q_(2n) from the formulas 138 to 143. The simultaneous equations are represented by the matrix of FIG. 11.

As described above, after deciding the arbitrary unknown constants p, q for the general solution regarding the voltage and current, the input/output voltage gain A_(v), the input impedance Z_(in) and the S parameter S for the transmission line 24 are obtained.

Firstly, for the input/output voltage gain A_(v), the voltage of each line segment θ is represented by the formula 114, and Z_(io) in the formula 134 is regarded as 0, whereby the input/output voltage gain A_(v) is obtained from the 144. $\begin{matrix} {A_{\upsilon} = \frac{\upsilon_{n}\left( {2\pi} \right)}{\upsilon_{1}(0)}} & \left( {{Formula}\quad 144} \right) \end{matrix}$

For the input impedance Z_(in), since Z_(1i) in the formula 132 is regarded as 0, the following formula 145 holds. If both sides of the formula 145 are divided by I₁(0), the following formula 146 results. Accordingly, the input impedance Z_(in) is obtained from the formula 147. $\begin{matrix} {{{V_{1}(0)} + {Z_{s} \cdot {I_{1}(0)}}} = V_{s}} & \left( {{Formula}\quad 145} \right) \\ {\frac{V_{s}}{I_{1}(0)} = {{Z_{s} + \frac{V_{1}(0)}{I_{1}(0)}}\quad = {Z_{s} + Z_{in}}}} & \left( {{Formula}\quad 146} \right) \\ {Z_{in} = \frac{V_{1}(0)}{I_{1}(0)}} & \left( {{Formula}\quad 147} \right) \end{matrix}$

Based on the output impedance Z, of the signal source V_(s) and the input impedance Z_(in) of the transmission line 24 with load Z_(t), the S₁₁ parameter is obtained from the following formula 148. Also, the S₂₁ parameter is obtained based on the S₁₁ parameter and the input/output voltage gain A_(v) from the following formula 149. $\begin{matrix} {S_{11} = \frac{Z_{in} - Z_{s}}{Z_{in} + Z_{s}}} & \left( {{Formula}\quad 148} \right) \\ {S_{21} = {\left( {1 + S_{11}} \right)A_{\upsilon}}} & \left( {{Formula}\quad 149} \right) \end{matrix}$

EXAMPLES

For a model including three lines 1 to 3 wired in parallel as shown in FIG. 12, a specific analysis example will be described. The lines 1 to 3 are magnetically coupled with each other, and capacitively coupled with each other. Also, in the case where they are capacitively coupled with the GND, the characteristics for high frequency three-phase AC current are analyzed. The specifications are as follows.

-   -   Length of line: 1.0 m     -   Signal sources V₁, V₂, V₃ (vector values) of lines 1, 2, 3:     -   V₁=5{cos(ωt)+j·sin(ωt)}     -   V₂=5 {cos(ωt+2/3π)+j·sin(ωt+2/3π)}     -   V₃=5 {cos(ωt−2/3π)+j·sin(ωt−2/3π)}     -   Output impedance of signal source: Z_(s1)=Z_(s2)=Z_(s3)=50 Ω     -   GND capacity: C₁=C₂=C₃=1000 pF     -   Capacity between lines: C₁₂=C₂₃=100 pF, C₃₁=50 pF     -   Distributed self-inductance of line: 25 μH/m     -   Connection coefficient: k=0.5     -   Line load: Z_(t1)=Z_(t2)=Z_(t3)=100 Ω

The unknown variables of this model (unknown variables p₁, q₁ of line 1, unknown variables p₂, q₂ of line 2, unknown variables p₃, q₃ of line 3) are decided. The boundary conditions for each of the lines 1 to 3 on the signal source are represented by the following formulas 150, 151 and 152. {dot over (V)} ₁ −i ₁(0)·Z _(s1)=ν₁(0)  (Formula 150) {dot over (V)} ₂ −i ₂(0)·Z _(s2)=ν₂(0)  (Formula 151) {dot over (V)} ₃ −i ₃(0)·Z_(s3)=ν₃(0)  (Formula 152)

These formulas are transformed into the following formulas 153, 154 and 155. ν₁(0)+i ₁(0)·Z _(s1) ={dot over (V)} ₁  (Formula 153) ν₂(0)+i ₂(0)·Z _(s2) ={dot over (V)} ₂  (Formula 154) ν₃(0)+i ₃(0)·Z _(s3) ={dot over (V)} ₃  (Formula 155)

The left sides of these formulas are rewritten with the general solution of the differential equation to obtain the following formulas 156, 157 and 158. $\begin{matrix} {{{\sum\limits_{i = 1}^{3}{u_{1i} \cdot p_{i}}} + {Z_{s1}{\sum\limits_{i = 1}^{3}{{\Psi\left( {i,1} \right)} \cdot q_{i}}}}} = {\overset{.}{V}}_{1}} & \left( {{Formula}\quad 156} \right) \\ {{{\sum\limits_{i = 1}^{3}{u_{2i} \cdot p_{i}}} + {Z_{s2}{\sum\limits_{i = 1}^{3}{{\Psi\left( {i,2} \right)} \cdot q_{i}}}}} = {\overset{.}{V}}_{2}} & \left( {{Formula}\quad 157} \right) \\ {{{\sum\limits_{i = 1}^{3}{u_{3i} \cdot p_{i}}} + {Z_{s3}{\sum\limits_{i = 1}^{3}{{\Psi\left( {i,3} \right)} \cdot q_{i}}}}} = {\overset{.}{V}}_{3}} & \left( {{Formula}\quad 158} \right) \end{matrix}$

On the other hand, the boundary conditions of each of the lines 1 to 3 on the load side are represented by the following formulas 159, 160 and 161. i ₁(2π)·Z _(l1)=ν₁(2π)  (Formula 159) i ₂(2π)·Z _(l2)=ν₂(2π)  (Formula 160) i ₃(2π)·Z _(l3)=ν₃(2π)  (Formula 161)

These formulas are transformed into the formulas 162, 163 and 164. ν₁(2π)−i ₁(2π)·Z _(l1)=0  (Formula 162) ν₂(2π)−i ₂(2π)·Z _(l2)=0  (Formula 163) ν₃(2π)−i ₃(2π)·Z _(l3)=0  (Formula 164)

The left sides of these formulas are rewritten with the general solution of the differential equation to obtain the following formulas 165, 166 and 167. $\begin{matrix} {{{\sum\limits_{i = 1}^{3}{\left\{ {{\cos\quad{h\left( {2\pi\quad\gamma_{i}} \right)}u_{1i}} - {Z_{l1}{\Psi\left( {i,1} \right)}\sin\quad{h\left( {2\pi\quad\gamma_{i}} \right)}}} \right\} p_{i}}} + {\sum\limits_{i = 1}^{3}{\left\{ {{\sin\quad{h\left( {2\pi\quad\gamma_{i}} \right)}u_{1i}} - {Z_{l1}{\Psi\left( {i,1} \right)}\cos\quad{h\left( {2\pi\quad\gamma_{i}} \right)}}} \right\} q_{i}}}} = 0} & \left( {{Formula}\quad 165} \right) \\ {{{\sum\limits_{i = 1}^{3}{\left\{ {{\cos\quad{h\left( {2\pi\quad\gamma_{i}} \right)}u_{2i}} - {Z_{l2}{\Psi\left( {i,2} \right)}\sin\quad{h\left( {2\pi\quad\gamma_{i}} \right)}}} \right\} p_{i}}} + {\sum\limits_{i = 1}^{3}{\left\{ {{\sin\quad{h\left( {2\pi\quad\gamma_{i}} \right)}u_{2i}} - {Z_{l2}{\Psi\left( {i,2} \right)}\cos\quad{h\left( {2\pi\quad\gamma_{i}} \right)}}} \right\} q_{i}}}} = 0} & \left( {{Formula}\quad 166} \right) \\ {{{\sum\limits_{i = 1}^{3}{\left\{ {{\cos\quad{h\left( {2\pi\quad\gamma_{i}} \right)}u_{3i}} - {Z_{l3}{\Psi\left( {i,3} \right)}\sin\quad{h\left( {2\pi\quad\gamma_{i}} \right)}}} \right\} p_{i}}} + {\sum\limits_{i = 1}^{3}{\left\{ {{\sin\quad{h\left( {2\pi\quad\gamma_{i}} \right)}u_{3i}} - {Z_{l3}{\Psi\left( {i,3} \right)}\cos\quad{h\left( {2\pi\quad\gamma_{i}} \right)}}} \right\} q_{i}}}} = 0} & \left( {{Formula}\quad 167} \right) \end{matrix}$

A total of six simultaneous equations composed the formulas 156, 157 and 158 and the formulas 165, 166 and 167 are solved to decide six unknown variables p₁, p₂, p₃, q₁, q₂, q₃.

Next, an output example will be described below. According to this invention, the distributions of the characteristic data of each of the lines 1 to 3, viz., voltage and current magnitude, phase, real element, imaginary element, and instantaneous value, are two-dimensionally displayed in different colors or densities corresponding to their values. As one example, FIG. 13A shows an output chart representing the voltage distribution of each of the lines 1 to 3 in densities with reference to a signal voltage (frequency 33 MHz), and FIG. 13B shows an output chart representing the phase distribution of each of the lines 1 to 3 in densities.

Also, in a three dimensional display in XYZ, it is possible that the distribution coordinates are allocated to two axes of X-Y, and the output magnitude is allocated to the remaining Z axis, for example. FIG. 14A shows a three dimensional display in which the voltage distribution of each of the lines 1 to 3 is represented in colors or densities, and FIG. 14B shows a three dimensional display in which the phase distribution of each of the lines 1 to 3 is represented in colors or densities.

Also, a graphical representation is allowed in which the frequency is taken along the transverse axis and the output value (or phase angle) is taken along the longitudinal axis. As one example, FIG. 15A shows the S₁₁ parameter characteristic of line 1, and FIG. 15B shows the S₂₁ parameter characteristic of line 1. In these figures, the bold line indicates the size of S parameter and the thin line indicates the phase angle of S parameter.

Also, a correlation graph between S₁₁ parameter and S₂₁ parameter may be displayed in a coordinate system in which the transverse axis is the size of S₁₁ parameter and the longitudinal axis is the size of S₂₁ parameter, as shown in FIG. 15C. Moreover, a vector diagram may be displayed in which the topology of S₁₁ parameter is allocated to the X vector component and the topology of S₂₁ parameter is allocated to the Y vector component, as shown in FIG. 15D. 

1. An apparatus for analyzing the distributed constant in a transmission line, comprising: an input section for inputting the distributed parameters for each line segment virtually and arbitrarily divided in said transmission line to be analyzed; an arithmetic operation section for calculating the characteristic data of said transmission line by performing a predetermined arithmetic operation based on said distributed parameters input from said input section; and an output section for outputting the characteristic data, wherein said distributed parameters input from said input section include a distributed self-inductance L_(i), a distributed resistance R_(i), and a distributed capacitance c_(i) and a distributed conductance g_(i) in relation to a reference potential plane in said transmission line, and a distributed mutual inductance M_(ij), a distributed capacitance C_(ij), and a distributed conductance G_(ij) between said transmission line and another transmission line, wherein said arithmetic operation section calculates the characteristic data in said transmission line by solving a differential equation regarding electromagnetic interaction of said line segment based on said predetermined distributed parameters.
 2. The apparatus for analyzing the distributed constant in the transmission line according to claim 1, wherein said arithmetic operation section comprises the secondary simultaneous differential equations as represented in the following formula 1 regarding electromagnetic interaction of said line segment with the voltage V(θ) of each line segment in said transmission line as an unknown variable (in the formula 1, V(θ) is a column vector of n×1 as represented in the following formula 2, Z is an impedance matrix created based on said distributed parameters as represented in the following formula 3, Y is an admittance matrix created based on said distributed parameters as represented in the following formula 4), in which said secondary simultaneous differential equations as represented in the formula 1 are solved by obtaining the eigenvalues γ_(i) ² (i=1 to n: n is the virtually divided number of line segments) of a complex coefficient matrix −ZY with said impedance matrix Z and said admittance matrix Y, thereby acquiring a general solution (p_(i), q_(i) are arbitrary constants, and u_(i) (i=1 to n) are eigenvectors corresponding to said eigenvalues γ_(i) ²) of voltage and current as represented in the following formula 5 and the following formula 6a or 6b at a limited position on said line segment. $\begin{matrix} {{\frac{\mathbb{d}^{2}}{\mathbb{d}\theta^{2}}{V(\theta)}} = {{- {ZY}} \cdot {V(\theta)}}} & \left( {{Formula}\quad 1} \right) \\ {{V(\theta)} = \begin{pmatrix} {\upsilon_{1}(\theta)} \\ {\upsilon_{2}(\theta)} \\ \vdots \\ {\upsilon_{n}(\theta)} \end{pmatrix}} & \left( {{Formula}\quad 2} \right) \\ {Z = \begin{pmatrix} \left( {R_{1} + {j\quad\omega\quad L_{1}}} \right) & {j\quad\omega\quad M_{12}} & \cdots & {j\quad\omega\quad M_{1n}} \\ {j\quad\omega\quad M_{21}} & \left( {R_{2} + {j\quad\omega\quad L_{2}}} \right) & \cdots & {j\quad\omega\quad M_{2n}} \\ \vdots & \vdots & ⋰ & \vdots \\ {j\quad\omega\quad M_{n1}} & {j\quad\omega\quad M_{n2}} & \cdots & \left( {R_{n} + {j\quad\omega\quad L_{n}}} \right) \end{pmatrix}} & \left( {{Formula}\quad 3} \right) \\ {{Y = \begin{pmatrix} \xi_{2} & \left( {G_{12} + {j\quad\omega\quad G_{12}}} \right) & \cdots & \left( {G_{1n} + {j\quad\omega\quad C_{1n}}} \right) \\ \left( {G_{21} + {j\quad\omega\quad C_{21}}} \right) & \xi_{2} & \cdots & \left( {G_{2n} + {j\quad\omega\quad C_{2n}}} \right) \\ \vdots & \vdots & ⋰ & \vdots \\ \left( {G_{n1} + {j\quad\omega\quad C_{n1}}} \right) & \left( {G_{n2} + {j\quad\omega\quad C_{n2}}} \right) & \cdots & \xi_{n} \end{pmatrix}}{{wherein},\quad{\xi_{j} = {- \left\{ {\left( {g_{j} + {j\quad\omega\quad c_{j}}} \right) + {\sum\limits_{{k = 1},{k \neq j}}^{n}\left( {G_{jk} + {j\quad\omega\quad C_{jk}}} \right)}} \right\}}}}} & \left( {{Formula}\quad 4} \right) \\ {{{V(\theta)} = {\sum\limits_{i = 2}^{n}{\left( {{p_{i}\quad\cos\quad{h\left( {\gamma_{i}\theta} \right)}} + {q_{i}\sin\quad{h\left( {\gamma_{i}\theta} \right)}}} \right)u_{i}}}}\quad} & \left( {{Formula}\quad 5} \right) \\ {{{I(\theta)} = {{Y{\int{{V(\theta)}{\mathbb{d}\theta}}}} = {Y{\sum\limits_{i = 1}^{n}{\frac{1}{\gamma_{i}}\left\{ {{p_{i}\quad\sin\quad{h\left( {\gamma_{i}\theta} \right)}} + {q_{i}\quad\cos\quad{h\left( {\gamma_{i}\theta} \right)}}} \right\} u_{i}}}}}}\quad{{or},}} & \left( {{Formula}\quad 6a} \right) \\ {{I(\theta)} = {{- Z^{1}}{\sum\limits_{i = 1}^{n}{\gamma_{i}\left\{ {{p_{1}\quad\sin\quad{h\left( {\gamma_{1}\theta} \right)}} + {q_{1}\cos\quad{{h1}\left( {\gamma_{1}\theta} \right)}}} \right\} d\quad\theta}}}} & \left( {{Formula}\quad 6b} \right) \end{matrix}$
 3. The apparatus for analyzing the distributed constant in the transmission line according to claim 2, wherein said arithmetic operation section further calculates arbitrary constants (unknown variables) p_(i), q_(i) (i=1 to n) included in said general solution of voltage and current as represented in said formula 5 and the formula 6a or 6b under the boundary conditions of said line segment that are input from said input section. $\begin{matrix} {{\frac{\mathbb{d}^{2}}{\mathbb{d}\theta^{2}}{I(\theta)}} = {{- {YZ}} \cdot {I(\theta)}}} & \left( {{Formula}\quad 7} \right) \\ {{I(\theta)} = {\sum\limits_{i = 1}^{n}{\left( {{p_{i}\cos\quad{h\left( {\gamma_{i}\theta} \right)}} + {q_{i}\sin\quad{h\left( {\gamma_{i}\theta} \right)}}} \right)u_{i}}}} & \left( {{Formula}\quad 8} \right) \\ {{V(\theta)} = {{{- Z}{\int{{I(\theta)}{\mathbb{d}\quad\theta}}}}\quad = {{- Z}{\sum\limits_{i = 1}^{n}{\frac{1}{\gamma_{i}}\left\{ {{p_{i}\quad\sin\quad{h\left( {\gamma_{i}\theta} \right)}} + {q_{i}{\cosh\left( {\gamma_{i}\theta} \right)}}} \right\} u_{i}}}}}} & \left( {{Formula}\quad 9a} \right) \\ {{or},} & \quad \\ {{V(\theta)} = {{- Z^{- 1}}{\sum\limits_{i = 1}^{n}{\gamma_{i}\left\{ {{p_{i}\quad\sin\quad{h\left( {\gamma_{i}\quad\theta} \right)}} + {q_{i}\quad\cos\quad{h\left( {\gamma_{i}\quad\theta} \right)}}} \right\} u_{i}}}}} & \left( {{Formula}\quad 9b} \right) \end{matrix}$
 4. The apparatus for analyzing the distributed constant in the transmission line according to claim 2, wherein said arithmetic operation section solves the secondary simultaneous differential equations as represented in the following formula 7 regarding electromagnetic interaction of said line segment with the current I(θ) of each line segment in said transmission line as an unknown variable, instead of said secondary simultaneous differential equations as represented in the formula 1, thereby acquiring a general solution of voltage and current as represented in the following formula 8 and the following formula 9a or 9b at a limited position on said line segment.
 5. The apparatus for analyzing the distributed constant in the transmission line according to claim 1, wherein said arithmetic operation section comprises the linearly simultaneous differential equations as represented in the following formula 10 regarding electromagnetic interaction of said line segment with the voltage V(θ) and current I(θ) of each line segment in said transmission line as unknown variables (in the formula 10, W(θ) is a matrix containing V(θ)_(n×1) and I(θ)_(n×1) as represented in the following formula 11, A is a matrix as represented in the following formula 12 containing said impedance matrix Z created based on said distributed parameters as represented in the formula 3 and said admittance matrix Y created based on said distributed parameters as represented in the formula 4), in which said linearly simultaneous differential equations as represented in the formula 10 are solved by obtaining 2n eigenvalues γ_(i) ²(i=1 to 2n) of said matrix A, thereby acquiring a general solution (α_(i)(i=1 to 2n) are arbitrary constants, and u_(i)(i=1 to 2n) are eigenvectors corresponding to said eigenvalues γ_(i)) of voltage and current as represented in the following formula
 13. $\begin{matrix} {{\frac{\mathbb{d}}{\mathbb{d}\theta}{W(\theta)}} = {A \cdot {W(\theta)}}} & \left( {{Formula}\quad 10} \right) \\ {{W(\theta)} = {\begin{pmatrix} {V(\theta)}_{n \times 1} \\ {I(\theta)}_{n \times 1} \end{pmatrix} = \begin{pmatrix} {\upsilon_{i}(\theta)} \\ {\upsilon_{2}(\theta)} \\ \vdots \\ {\upsilon_{n}(\theta)} \\ {i_{1}(\theta)} \\ {i_{2}(\theta)} \\ \vdots \\ {i_{n}(\theta)} \end{pmatrix}}} & \left( {{Formula}\quad 11} \right) \\ {A = \begin{pmatrix} 0_{n \times n} & {- Z} \\ Y & 0_{n \times n} \end{pmatrix}_{2n \times 2n}} & \left( {{Formula}\quad 12} \right) \\ {{W(\theta)} = {\sum\limits_{i = 1}^{2n}{\alpha_{i}{\mathbb{e}}^{\gamma\quad i\quad\theta}u_{i}}}} & \left( {{Formula}\quad 13} \right) \end{matrix}$
 6. The apparatus for analyzing the distributed constant in the transmission line according to claim 5, wherein said arithmetic operation section sets up 2n simultaneous equations under the boundary conditions of said line segment input from said input section for each of the j-th (j=1 to n) components of a voltage vector as represented in the following formula 14 and a current vector as represented in the following formula 15 in the transmission line, which are derived from the general solution as represented in the formula 13, thereby acquiring 2n arbitrary constants (unknown variables) α_(i). $\begin{matrix} {{v_{j}(\theta)} = {\sum\limits_{i = 1}^{2n}{\alpha_{i}{\mathbb{e}}^{\gamma\quad i\quad\theta}u_{ji}}}} & \left( {{Formula}\quad 14} \right) \\ {{i_{j}(\theta)} = {\sum\limits_{i = 1}^{2n}{\alpha_{i}{\mathbb{e}}^{\gamma\quad i\quad\theta}u_{n + {ji}}}}} & \left( {{Formula}\quad 15} \right) \end{matrix}$
 7. The apparatus for analyzing the distributed constant in the transmission line according to claim 3, wherein said arithmetic operation section calculates an input/output voltage gain A_(v) from an input voltage and a terminating voltage of said signal line after obtaining the arbitrary constants of the general solution regarding the voltage and current.
 8. The apparatus for analyzing the distributed constant in the transmission line according to claim 3, wherein said arithmetic operation section calculates an input impedance Z_(in) of said signal line based on a signal source voltage V_(s) and a signal source output impedance Z_(s) included in the boundary conditions given from said input section after obtaining the arbitrary constants of the general solution regarding the voltage and current.
 9. The apparatus for analyzing the distributed constant in the transmission line according to claim 3, wherein said arithmetic operation section calculates an input/output voltage gain A_(v) from an input voltage and a terminating voltage of said signal line after obtaining the arbitrary constants of the general solution regarding the voltage and current, and calculates an input impedance Z_(in) of said signal line based on a signal source voltage V_(s) and a signal source output impedance Z_(s) included in the boundary conditions given from said input section, thereby acquiring an S parameter S₁₁ from said signal source output impedance Z_(s) and said input impedance Z_(in), and an S parameter S₂₁ from said S parameter S₁₁ and said input/output voltage gain A_(v).
 10. The apparatus for analyzing the distributed constant in the transmission line according to claim 1, wherein said signal line is involved in the equipment such as a connector, a sensor or a filter to allow the distributed constant characteristics of said equipment itself to be analyzed.
 11. The apparatus for analyzing the distributed constant in the transmission line according to claim 2, wherein said arithmetic operation section solves the secondary simultaneous differential equations as represented in the formula 1 with the voltage V(θ) of each line segment as an unknown variable under the boundary conditions regarding the voltage of said each line segment given from said input section in accordance with a numerical solution algorithm of the differential equation.
 12. The apparatus for analyzing the distributed constant in the transmission line according to claim 4, wherein said arithmetic operation section solves the secondary simultaneous differential equations as represented in the formula 7 with the current I(θ) of each line segment as an unknown variable under the boundary conditions regarding the current of said line segment given from said input section in accordance with a numerical solution algorithm of the differential equation.
 13. The apparatus for analyzing the distributed constant in the transmission line according to claim 5, wherein said arithmetic operation section solves the linearly simultaneous differential equations as represented in the formula 10 with the voltage V(θ) and current I(θ) of each line segment as unknown variables under the boundary conditions regarding the voltage and current of said line segment given from said input section in accordance with a numerical solution algorithm of the differential equation.
 14. The apparatus for analyzing the distributed constant in the transmission line according to claim 11, wherein said numerical solution algorithm is any one of the Runge-Kutta methods, Shooting methods (Gun methods), and Relaxation methods.
 15. A program for use with an apparatus for analyzing the distributed constant in a transmission line, wherein said arithmetic operation section of claim 2, which consists of a computer, sets up the secondary simultaneous differential equations as represented in the following formula 1 regarding electromagnetic interaction of said line segment with the voltage V(θ) of each line segment virtually and arbitrarily divided in said transmission line to be analyzed as an unknown variable (in the formula 1, V(θ) is a column vector of n×1 as represented in the following formula 2, Z is an impedance matrix created based on said distributed parameters as represented in the following formula 3, Y is an admittance matrix created based on said distributed parameters as represented in the following formula 4), and solves said secondary simultaneous differential equations as represented in the formula 1 by obtaining the eigenvalues γ_(i) ² (i=1 to n:n is the virtually divided number of line segments) of a complex coefficient matrix −ZY with said impedance matrix Z and said admittance matrix Y from a distributed self-inductance L_(i), a distributed resistance R_(i), and a distributed capacitance c_(i) and a distributed conductance g_(i) in relation to a reference potential plane in said transmission line, and a distributed mutual inductance M_(ij), a distributed capacitance C_(ij), and a distributed conductance G_(ij) between said transmission line and another transmission line, which are input from said input section, thereby acquiring a general solution (p_(i), q_(i) are arbitrary constants, and u; (i=1 to n) are eigenvectors corresponding to said eigenvalues γ_(i) ²) of voltage and current as represented in the following formula 5 and the following formula 6a or 6b at a limited position on said line segment. $\begin{matrix} {{\frac{\mathbb{d}^{2}}{\mathbb{d}\theta^{2}}{V(\theta)}} = {{- {ZY}} \cdot {V(\theta)}}} & \left( {{Formula}\quad 1} \right) \\ {{V(\theta)} = \begin{pmatrix} {\upsilon_{1}(\theta)} \\ {\upsilon_{2}(\theta)} \\ \vdots \\ {\upsilon_{n}(\theta)} \end{pmatrix}} & \left( {{Formula}\quad 2} \right) \\ {Z = \begin{pmatrix} \left( {R_{1} + {j\quad\omega\quad L_{1}}} \right) & {j\quad\omega\quad M_{12}} & \cdots & {j\quad\omega\quad M_{1n}} \\ {j\quad\omega\quad M_{21}} & \left( {R_{2} + {j\quad\omega\quad L_{2}}} \right) & \cdots & {j\quad\omega\quad M_{2n}} \\ \vdots & \vdots & ⋰ & \vdots \\ {j\quad\omega\quad M_{n1}} & {j\quad\omega\quad M_{n2}} & \cdots & \left( {R_{n} + {j\quad\omega\quad L_{n}}} \right) \end{pmatrix}} & \left( {{Formula}\quad 3} \right) \\ {{Y = \begin{pmatrix} \xi_{1} & \left( {G_{12} + {j\quad\omega\quad C_{12}}} \right) & \cdots & \left( {G_{1n} + {j\quad\omega\quad C_{1n}}} \right) \\ \left( {G_{21} + {j\quad\omega\quad C_{21}}} \right) & \xi_{2} & \cdots & \left( {G_{2n} + {j\quad\omega\quad C_{2n}}} \right) \\ \vdots & \vdots & ⋰ & \vdots \\ \left( {G_{n1} + {j\quad\omega\quad C_{n1}}} \right) & \left( {G_{n2} + {j\quad\omega\quad C_{n2}}} \right) & \cdots & \xi_{n} \end{pmatrix}}{{wherein},\quad{\xi_{j} = {- \left\{ {\left( {g_{j} + {j\quad\omega\quad c_{j}}} \right) + {\sum\limits_{{k = 1},{k \neq j}}^{n}\left( {G_{jk} + {j\quad\omega\quad C_{jk}}} \right)}} \right\}}}}} & \left( {{Formula}\quad 4} \right) \\ {{{V(\theta)} = {\sum\limits_{i = 1}^{n}{\left( {{p_{i}\quad\cos\quad{h\left( {\gamma_{i}\theta} \right)}} + {q_{i}\sin\quad{h\left( {\gamma_{i}\theta} \right)}}} \right)u_{i}}}}\quad} & \left( {{Formula}\quad 5} \right) \\ {{{I(\theta)} = {{Y{\int{{V(\theta)}{\mathbb{d}\theta}}}} = {Y{\sum\limits_{i = 1}^{n}{\frac{1}{\gamma_{i}}\left\{ {{p_{i}\quad\sin\quad{h\left( {\gamma_{i}\theta} \right)}} + {q_{i}\quad\cos\quad{h\left( {\gamma_{i}\theta} \right)}}} \right\} u_{i}}}}}}\quad{{or},}} & \left( {{Formula}\quad 6a} \right) \\ {{I(\theta)} = {{- Z^{- 1}}{\sum\limits_{i = 1}^{n}{\gamma_{i}\left\{ {{p_{1}\quad\sin\quad{h\left( {\gamma_{i}\theta} \right)}} + {q_{1}\cos\quad{h\left( {\gamma_{1}\theta} \right)}}} \right\} d\quad\theta}}}} & \left( {{Formula}\quad 6b} \right) \end{matrix}$
 16. The program for use with the apparatus for analyzing the distributed constant in the transmission line according to claim 15, wherein said computer (arithmetic operation section) further calculates arbitrary constants (unknown variables) p_(i), q_(i) (i=1 to n) included in the general solution of voltage and current as represented in the formula 5 and the formula 6a or 6b under the boundary conditions of said line segment that are input from said input section.
 17. The program for use with the apparatus for analyzing the distributed constant in the transmission line according to claim 15, wherein said computer (arithmetic operation section) sets up the secondary simultaneous differential equations as represented in the following formula 7 regarding electromagnetic interaction of said line segment with the current I(θ) of each line segment in said transmission line as an unknown variable, instead of said secondary simultaneous differential equations as represented in the formula 1; and solves said secondary simultaneous differential equations as represented in the formula 7 by obtaining the eigenvalues γ_(i) ²(i=1 to n: n is the virtually divided number of line segments) of a complex coefficient matrix −ZY with said impedance matrix Z and said admittance matrix Y, thereby acquiring a general solution of voltage and current as represented in the following formula 8 and the following formula 9a or 9b at a limited position on said line segment. $\begin{matrix} {{\frac{\mathbb{d}^{2}}{\mathbb{d}\theta^{2}}{I(\theta)}} = {{- {YZ}} \cdot {I(\theta)}}} & \left( {{Formula}\quad 7} \right) \\ {{I(\theta)} = {\sum\limits_{i = 1}^{n}{\left( {{p_{i}\cos\quad{h\left( {\gamma_{i}\theta} \right)}} + {q_{i}\sin\quad{h\left( {\gamma_{i}\theta} \right)}}} \right)u_{i}}}} & \left( {{Formula}\quad 8} \right) \\ {{V(\theta)} = {{{- Z}{\int{{I(\theta)}{\mathbb{d}\quad\theta}}}}\quad = {{- Z}{\sum\limits_{i = 1}^{n}{\frac{1}{\gamma_{i}}\left\{ {{p_{i}\quad\sin\quad{h\left( {\gamma_{i}\theta} \right)}} + {q_{i}{\cosh\left( {\gamma_{i}\theta} \right)}}} \right\} u_{i}}}}}} & \left( {{Formula}\quad 9a} \right) \\ {{or},} & \quad \\ {{V(\theta)} = {{- Z^{- 1}}{\sum\limits_{i = 1}^{n}{\gamma_{i}\left\{ {{p_{i}\quad\sin\quad{h\left( {\gamma_{i}\quad\theta} \right)}} + {q_{i}\quad\cos\quad{h\left( {\gamma_{i}\quad\theta} \right)}}} \right\} u_{i}}}}} & \left( {{Formula}\quad 9b} \right) \end{matrix}$
 18. A program for use with the apparatus for analyzing the distributed constant in the transmission line, wherein said arithmetic operation section of claim 1, which consists of a computer, sets up the linearly simultaneous differential equations as represented in the following formula 10 regarding electromagnetic interaction of said line segment with the voltage V(θ) and current I(θ) of each line segment virtually and arbitrarily divided in said transmission line to be analyzed as unknown variables (in the formula 10, W(θ) is a matrix containing V(θ)_(n×1) and I(θ)_(n×1) as represented in the following formula 11, A is a matrix as represented in the following formula 12 containing said impedance matrix Z created based on said distributed parameters as represented in the formula 3 and said admittance matrix Y created based on said distributed parameters as represented in the formula 4), and solves said linearly simultaneous differential equations as represented in the formula 10 by obtaining 2n eigenvalues γ_(i)(i=1 to 2n) of said matrix A from a distributed self-inductance L_(i), a distributed resistance R_(i), and a distributed capacitance ci and a distributed conductance g_(i) in relation to a reference potential plane in said transmission line, and a distributed mutual inductance M_(ij), a distributed capacitance C_(ij), and a distributed conductance G_(ij) between said transmission line and another transmission line, which are input from said input section, thereby acquiring a general solution (α_(i)(i=1 to 2n) are arbitrary constants, and u_(i)(i=1 to 2n) are eigenvectors corresponding to said eigenvalues γ_(i)) of voltage and current as represented in the following formula
 13. $\begin{matrix} {{\frac{\mathbb{d}}{\mathbb{d}\theta}{W(\theta)}} = {A \cdot {W(\theta)}}} & \left( {{Formula}\quad 10} \right) \\ {{W(\theta)} = {\begin{pmatrix} {V(\theta)}_{n \times 1} \\ {I(\theta)}_{n \times 1} \end{pmatrix} = \begin{pmatrix} {\upsilon_{i}(\theta)} \\ {\upsilon_{2}(\theta)} \\ \vdots \\ {\upsilon_{n}(\theta)} \\ {i_{1}(\theta)} \\ {i_{2}(\theta)} \\ \vdots \\ {i_{n}(\theta)} \end{pmatrix}}} & \left( {{Formula}\quad 11} \right) \\ {A = \begin{pmatrix} 0_{n \times n} & {- Z} \\ Y & 0_{n \times n} \end{pmatrix}_{2n \times 2n}} & \left( {{Formula}\quad 12} \right) \\ {{W(\theta)} = {\sum\limits_{i = 1}^{2n}{\alpha_{i}{\mathbb{e}}^{\gamma\quad i\quad\theta}u_{i}}}} & \left( {{Formula}\quad 13} \right) \end{matrix}$
 19. The program for use with the apparatus for analyzing the distributed constant in the transmission line according to claim 18, wherein said computer (arithmetic operation section) sets up 2n simultaneous equations under the boundary conditions of said line segment input from said input section for each of the j-th (j=1 to n) components of a voltage vector as represented in the following formula 14 and a current vector as represented in the following formula 15 in the transmission line, which are derived from the general as represented in the formula 13, thereby acquiring 2n arbitrary constants (unknown variables) α_(i). $\begin{matrix} {{v_{j}(\theta)} = {\sum\limits_{i = 1}^{2n}{\alpha_{i}{\mathbb{e}}^{\gamma\quad i\quad\theta}u_{ji}}}} & \left( {{Formula}\quad 14} \right) \\ {{i_{j}(\theta)} = {\sum\limits_{i = 1}^{2n}{\alpha_{i}{\mathbb{e}}^{\gamma\quad i\quad\theta}u_{n + {ji}}}}} & \left( {{Formula}\quad 15} \right) \end{matrix}$ 